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In this communication, we focus on the possibility to test GR with radioscience experiments. 
We present a new software that in a first step simulates the Range/Doppler signals directly 
from the space time metric (thus in GR and in alternative theories of gravity). In a second 
step, a least-squares fit of the involved parameters is performed in GR. This software allows one 
to get the order of magnitude and the signature of the modifications induced by an alternative 
theory of gravity on radioscience signals. As examples, we present some simulations for the 
Gassini mission in Post-Einsteinian gravity and with the MOND External Field Effect. 



1 Introduction 

There is still a great interest in testing General Relativity (theoretical motivations such as quan- 
tification of gravity, unification with other interactions. . . ). Within the solar system, the grav- 
itational observations are always related with radioscience measurements (Range and Doppler) 
or with angular measurements (position of body in the sky, VLBI). In this communication, we 
present a new tool that performs simulations of radioscience experiments in General Relativity 
and in alternative metric theories of gravitation. The developed software simulates Range and 
Doppler signals directly from the space time metric (and from the initial conditions of the bod- 
ies considered). The aim of this tool is to provide orders of magnitude and signatures for the 
variations of the Range/Doppler induced by alternative theories of gravity. In order to compare 
signals in difi'erent theories, a least-squares fit of the different parameters involved in the problem 
(mainly the initial conditions) is performed. In this communication, we present the principles 
of the software and results of simulations of the Cassini probe missio n du ring its cruise from 
Jupiter to Saturn with Post-Einsteinian modifications of gravity (PEG)^^] or with an external 
field effect due to a MONDian modification of gravity^. 



2 Description of the software 

Our software directly generates Doppler/Range signal from the metric considered in a fully- 
relativistic way. This includes the integration of the equations of motion in a given coordinate 
time, the computation of the time transfer of light and the clock behavior. 



The equations of motion are derived from the metric and the geodesic equations integrated 
with respect to coordinate time^. 

The connection between the coordinate and proper time of clocks is obtained by integrating 
the equation of proper time^ 

^ = y'soo - '^goiVi - gijViVj (1) 

where g^u is the space-time metric, Vi is the velocity of the clock and r its proper time. 

Finally, the time transfer is also determined directly from the metric using Synge World's 
function formalism^. Within this formalism, one does not need to integrate the photon tra- 
jectory in order to get the time transfer or the frequency shift (in the linear approximation). 
Instead, those quantities are expressed as integrals of functions defined from the metric (and 
their derivatives) along the photon Minkowski path. For example the coordinate propagation 
time can be expressed as ^ 



with 



T{xlite),xl{tr),tr) = ^ + ^ ( f{z^{f^))dfl (2) 

c c Jo 



f = -hoo-2Ni,ho,-Ni,Ni>,hij, h^,y = g^,y-ri^y, AT*^ = (3) 



xl,te are the position and time of the emitter (computed iteratively assuming flat space time), 
are the position and time of the receptor and i?er = ||a^e(*e) — a^rC^r)]]- The integral in ([2]) 
is performed over a Minkowski path between emitter and receptor (z"(/i) is a straight line). A 
similar expression is used for the frequency shift. 

In order to investigate the observable signatures of an alternative theory of gravity in the 
Range and Doppler data we perform a least-squares fit in GR on the different parameters (initial 
conditions and masses of the bodies) and search for identifiable signatures in the residuals. The 
fit of initial conditions is necessary in order to avoid effects due to the choice of coordinates. As 
a matter of fact, this fit is always done in practice, which means that a reasonable analysis of 
the signal produced by an alternative theory has to be done after the fit. 



3 Simulations of Cassini mission 



As an example, we simulate the two-way range and Doppler signals to the Cassini spacecraft 
from June 2002 during 3 years (when the probe was between Jupiter and Saturn). To simplify 
the situation we only consider the Sun, the Earth and Cassini spacecraft. 



3.1 Post-Einsteinian Gravity (PEG) 

The first alternative metric theory considered is Post-Einsteinian Gravity (PEG)II^. From a 
phenomenological point of view this theory consists in including two potentials $Ar(r) and $p(r) 
to the metric. Here we concentrate on the sector <I>p(r). We consider a series expansion. That 
is to say we suppose the spatial part of the metric to be modified as follows 

9ij = [9ij]GR - (xir + X2r^ + '^7^ j (4) 

where r is a radial isotropic coordinate, M is the sun mass, c the velocity of light and G the 
gravitational constant. The parameter 6'j is related to the post-newtonian parameter ^7 = 7 — 1. 

Different simulations were performed with different values of the three PEG parameters. 
For example. Figure [J represents the Range and Doppler differences between a simulation in a 



theory with ^7 = 7 — 1 = 10~^ and m GR. The three peaks occur during solar conjunctions. 
The signal due to the conjunction is not absorbed at all by the fit of the initial conditions which 
nevertheless absorbs a large modulations in the range signal. 
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Figure 1: Representation of the Range (on the left) and Doppler (on the right) signals due to an alternative theory 
with 7 — 1 = 10^''. The blue line is the difference between a simulation in the alternative theory and a simulation 
in GR (with the same parameters) . The green line is the residuals obtained after analyzing the simulated data in 
GR (which means after the fit of the different parameters). 



To summarize, Figure [5] represents the maximal difference between the Doppler generated 
in PEG theory and the Doppler generated in GR for different PEG theories (characterized by 
their values of xij Xi ^"^^ ^l)- If we request the residuals to be smaller than Cassini Doppler 
accuracy (roughly 10~^^), we get boundary values for the three parameters: xi < 5 10~'^^m~^, 
X2 < 2 10~^'^m~^ and 7 — 1 < 3 10~^ (which is very similar to the real estimation!^. 




Figure 2: Representation of the maximal Doppler signal due to PEG theory (parametrized by three parameters 
XXiXiil — 1) for the Cassini mission between Jupiter and Saturn. The blue lines represent the maximal Doppler 
difference between a simulation in the alternative theory and a simulation in GR (with the same parameters). 
The green lines represent the maximal residuals obtained after analyzing the simulated data's in GR (i.e. after 
the fit of the parameters). The read lines represent the assumed Cassini accuracy. 



3.2 MOND External Field Effect (EFE) 

Another alternative theory considered is the External Field Effect produced by a MOND the- 
ory^. In this framework, the dominant effect is modeled by a quadrupolar contribution to the 
Newtonian potential U = + ^^x'^x^ i^^i^j ~ where is a unitary vector pointing to- 

wards the galactic center and 2.1 10^^''s^^ ^ Q2 1^ 4.1 lO^^^s^^ is the value of the quadrupole 
moment whose value depends on the MOND function. 



Figure [3] represents the effect of the EFE on the Range and Doppler signals from Cassini. It 
can be seen that the signal are just below the Cassini accuracy (10~^^ in Doppler). Therefore, 
the Cassini arc considered here is not sufficient to provide a good test of MOND External Field 
Effect. 
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Figure 3: Representation of the Range (on the left) and Doppler (on the right) signals due to the MOND EFE 
(Q2 = 4.1 10~^^s~^). The blue line is the difference between a simulation with the EFE and a simulation in GR 
(with the same parameters). The green line shows the residuals obtained after analyzing the simulated data in 
GR (which means after the fit of the different parameters). 



4 Conclusion 

In this communication, we have presented a new tool that performs Range/Doppler simulations 
in metric theories of gravity. With this software, it is easy to get the order of magnitude and the 
signature of the modifications induced by alternative theories of gravity on radioscience signals. 
As an example, we have presented some simulations for the Cassini mission in Post-Einsteinian 
Gravity and we have derived boundary values for some PEG parameters. We have also presented 
simulations including the MOND External Field Effect and we have shown that this effect is too 
small to be detected during the Cassini cruise between Jupiter and Saturn. 

In the future, further simulations can be done for other theories and other (future and past) 
space missions. 
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